Simultaneous inversion of mantle properties and initial conditions using an adjoint of mantle convection

Through the assimilation of present-day mantle seismic structure, adjoint methods can be used to constrain the structure of the mantle at earlier times, i.e., mantle initial conditions. However, the application to geophysical problems is restricted through both the high computational expense from repeated iteration between forward and adjoint models and the need to know mantle properties (such as viscosity and the absolute magnitude of temperature or density) a priori. We propose that an optimal first guess to the initial condition can be obtained through a simple backward integration (SBI) of the governing equations, thus lessening the computational expense. Given a model with known mantle properties, we show that a solution based on an SBI-generated first guess has smaller residuals than arbitrary guesses. Mantle viscosity and the effective Rayleigh number are crucial for mantle convection models, neither of which is exactly known. We place additional constraints on these basic mantle properties when the convection-induced dynamic topography on Earth's surface is considered within an adjoint inverse method. Besides assimilating present-day seismic structure as a constraint, we use dynamic topography and its rate of change in an inverse method that allows simultaneous inversion of the absolute upper and lower mantle viscosities, scaling between seismic velocity and thermal anomalies, and initial condition. The theory is derived from the governing equations of mantle convection and validated by synthetic experiments for both one-layer viscosity and two-layer viscosity regionally bounded spherical shells. For the one-layer model, at any instant of time, the magnitude of dynamic topography is controlled by the temperature scaling while the rate of change of topography is controlled by the absolute value of viscosity. For the two-layer case, the rate of change of topography constrains upper mantle viscosity while the magnitude of dynamic topography determines the temperature scaling (lower mantle viscosity) when upper-mantle (lower-mantle) density anomaly dominates the flow field; this two-stage scheme minimizes the tradeoff between temperature and lower mantle viscosity. For both cases, we show that the theory can constrain mantle properties with errors arising through the adjoint recovery of the initial condition; for the two-layer model, this error is manifest as a tradeoff between the temperature scaling and lower mantle viscosity

Adjoint models of mantle convection with seismic, plate motion, and stratigraphic constraints: North America since the Late Cretaceous

We apply adjoint models of mantle convection to North America since the Late Cretaceous. The present-day mantle structure is constrained by seismic tomography and the time-dependent evolution by plate motions and stratigraphic data (paleoshorelines, borehole tectonic subsidence, and sediment isopachs). We infer values of average upper and lower mantle viscosities, provide a synthesis of North American vertical motions (relative sea level) from the Late Cretaceous to the present, and reconstruct the geometry of the Farallon slab back to the Late Cretaceous. In order to fit Late Cretaceous marine inundation and borehole subsidence, the adjoint model requires a viscosity ratio across 660 km discontinuity of 15:1 (reference viscosity of 10^(21) Pa s), which is consistent with values previously inferred by postglacial rebound studies. The dynamic topography associated with subduction of the Farallon slab is localized in western North America over Late Cretaceous, representing the primary factor controlling the widespread flooding. The east coast of the United States is not stable; rather, it has been experiencing continuous dynamic subsidence over the Cenozoic, coincident with an overall eustatic fall, explaining a discrepancy between sea level derived from the New Jersey coastal plain and global curves. The east coast subsidence further constrains the mantle viscosity structure and requires an uppermost mantle viscosity of 10^(20) Pa s. Imposed constraints require that the Farallon slab was flat lying during Late Cretaceous, with an extensive zone of shallow dipping Farallon subduction extending beyond the flat-lying slab farther east and north by up to 1000 km than previously suggested